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ABSTRACT 

A novel model for collimation and transport of electron-positron-ion jets is presented. 
Analytical results show that the filamentary structures can be sustained by self-induced toroidal 
magnetic fields permeating through the filaments, whose widths significantly expand in the 
pair-dominant regimes. The magnetic field strength reflects a characteristic of equipartition 
of excess kinetic energy of the jets. It is also shown that growth of the hoselike instability 
is strongly suppressed. Essential features derived from this model are consistent with recent 
results observed by using the very long baseline telescopes. 

Subject headings: galaxies: jets — magnetic fields — methods: analytical — plasmas 



1. INTRODUCTION 

Now, there still exist some challenging problems unsolved on astrophysical jets: these include the 
acceleration (Blandford & Payne 1982), collimation, and stability (Benford 1978; Appl & Camenzind 
1992). For the past decades, substantial studies have been devoted to investigating the directional mass 
acceleration based on the magnetohydrodynamics (MHD; Shibata & Uchida 1986, 1990; Cao & Spruit 
1994; Bell 1994; Matsumoto et al. 1996; Kudoh & Shibata 1997a, 1997b) and radiative (Tajima & Fukue 
1996, 1998; Fukue, Tojyo, & Hirai 2001) mechanism around the central engine of the accretion disk. Koide, 
Nishikawa, & Mutel (1996) have performed the numerical simulation of a relativistic MHD (RMHD) jet 
injected into a poloidal magnetic field. The three-dimensional RMHD simulations are in progress worldwide 
(Nishikawa et al. 1997; Aloy et al. 1999). In addition, Carilli & Barthel (1996) have proposed the ballistic, 
pressure-matched jets whose internal pressure was assumed to be comparable to or even less than the 
external one. 

On the other hand, one can find observational evidence that the energy density of jets is larger than that 
of the cocoon or intergalactic medium (e.g., 4C 32.69; Pottash & Wardle 1980). Recent rotation measures 
using very long baseline interferometry (VLBI) also reveal that the direction of magnetic field vectors is 
likely to be transverse to the jet axis (e.g., 1803+784; Gabuzda 1999). In some cases, the transverse fields 
are so smooth that one can barely discriminate the knots, and it seems difficult to wholly understand such 
polarization properties by invoking a trail of oblique compressional shocks. A puzzling question is, therefore, 
how to demonstrate the orientation of magnetic fields as well as the related se//-collimation mechanism, 
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which must sustain the large-scale structure as observed. Nevertheless, to the best of our knowledge, there 
is no publication that proposes such a coherent scenario in the context of astrophysical jets, especially 
active galactic nucleus (AGN) jets, which extend up to megaparsec scales (Mpc ~ 3 x 10 24 cm) with narrow 
opening angles {l<j> <; 1°) in some objects (e.g., NGC 6251: Bridle & Perley 1984; Cyg A: Carilli et al. 
1998; M87, PKS 0637-752: Hirabayashi 2001). 

We address in this Letter that the toroidal (transverse) magnetic field induced by electron-positron 
flow is a favored candidate for the collimation force of the AGN jets. We find that the screening effects 
of electron-positron-ion plasmas play a significant role in the self-organization of the filament envelopes, 
whose radii exceedingly stretch in the positron-rich regimes. The screening by the pairs and the diffuse 
envelope prevent jets from the snakelike distortion. A scaling law of the maximum possible magnetic fields 
is displayed as a matter of convenience. We expect that the present model is also applicable to Galactic 
jets, except for neutral flows from protostars, and so on (Tajima & Shibata 1997). 

The key issue discussed here is the recurrence of the "plasma universe" model proposed by Alfven 
(1981), which is relevant to cosmic ray generation and transport involving pinches (Trubnikov 1991). So far, 
the present plasma configuration itself has been considered to be notoriously unstable, but we attempt to 
rewrite that scenario. The new points, going beyond the pioneering work by Benford (1978), are mainly (1) 
the pairs are introduced, (2) the scaling of jet radius and field strength is presented, and (3) detailed kinetic 
theory is applied to the instability analysis. We also mention that the underlying physics is ubiquitous. 
Regarding recent laboratory experiments, the high power fusion lasers up to a petawatt (10 15 W) reproduce 
the collimated relativistic electron jets (Key et al. 1998), several tens MeV ions (Snavely et al. 2000), 
and even the Bethe-Heitler pair productions triggered by the bremsstrahlung 7-ray photons (Cowan et al. 
1999). These experiments are now opening door to the high-field laboratory astrophysics. 

As ordinarily expected, in fully a relativistic regime of T > 10 10 K, ejecta could consist of electrons, 
positrons, and a small portion of ions, i.e., n e -^>n e+ ^$>rii, where n e -, n e+ , and rij are the number densities 
of electrons, positrons, and ions, respectively. Recently, the ASCA satellite has also detected X-ray emissions 
of various ion species, such as Fe, Ni, Mg, Si, S, Nc, and Ar, from the SS 433 jets (Kotani et al. 1996). 
In optically thick regions, the electron-positron photoplasmas may establish the quasi- Wien equilibria 
(Iwamoto & Takahara 2002). Anyhow, around central engines, fast electron-positron flows relative to ion 
motions can be easily organized owing to the difference of their inertia. Namely, the electron-positron clouds 
are primarily accelerated by, e.g., the Lorentz force: | =Fe(w eT /c) x B\/m eT » \{Z*)e(vi on /c) x B\/((A)m p ), 
where m p /m eT — 1836, and (Z*) and (A) are the averages of charge state and mass number of the 
multispecies ions, respectively. 

The most important point is that the stream along the jet axis z with the excess electrons 
n*_ = n e _ — n e+ = {Z*)ni generates the toroidal "self -magnetic field Bg, which participates in 
self-pinching the electron-positron gas and in assembling the ions radially inwards on the hydrodynamical 
timcscale (Honda, Meyer-ter-Vehn, & Pukhov 2000a, 2000b). Since the plasma holds a quite high 
conductivity, a closed current system including a return current inside and/or outside the jets is 
self-organized without delay. It is noted is that the toroidal magnetic fields act as a defocusing force for 
the forward-running beam positrons and backward return electrons, whereas such deflections tend to be 
immediately restored by the electrostatic fields due to microscopic charge separations, to avoid unphysical 
charge-up. Therefore, one can treat the electron-positron flow as the negatively charged fluid, which 
partially compensates for the positive ion background, as shown below. 
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2. THE ELECTRON-POSITRON FLUID EQUATIONS 



We start with a relativistic electron-positron beam-plasma equilibrium, invoking that the nonrelativistic 
equation of state of ideal gas p — nT is valid for the relativistic bulk motion (Rindlcr 1982). Taking 
the return components into account, radial force balance on the electron-positron fluid elements can be 
expressed as dpj T /dr = q T rij T (E r + Vj T xBg/c), where j = b,r indicate the beam and return components, 
respectively. For electrons and positrons, we set q- = — |e and q + = e|, respectively. The photon 
pressure is omitted, and the optically thin jet is sufficiently distant from the central engine. All physical 
quantities are in the ion rest (jet) frame throughout this Letter; this frame is chosen even in the case of 
n e _ ss n e+ > 2 x 10 3 (A)n i . The force balance equations can be then cast to 



d_ 

dr 



( / 6 n e _(r)T 6 _ 
(1 - / 6 )n e -(r)T r _ 
f b n e+ (r)T b+ 
V (1 - fb)n e +(r)T r+ 



( -ef b n e -(r) 
-e(l - f b )n e -(r) 
efbn e +(r) 
V e(l - f b )n e+ (r) 



e(3 b -f b n e -(r) \ 
e/3 r _(l - f b )n e -(r) 

-ef3 b+ f b n e+ (r) 
-e(3 r+ (l - f b )n e+ {r) ) 



E r (r) 



(1) 



/ x 47T 

Er(r) = — 
r 



Vg T / dr'r'n jT (r') + (Z*)e f dr'r'n^r') 



(2) 



B s e (r) = - E «t /' dr'r'P jT n jT (r'), (3) 

where n bT (r) = f b n eT (r) and n rT (r) = (1 — f b )n eT (r) reflect the portions of the forward beam and 
backward return flows, respectively. Furthermore, we introduce the ratio of positron/electron densities, 



< f P 



/n e 



const < 1 and the fractional charge neutrality of ion/electron charge densities, 



< f c = (Z*)rii/(n e - — n e+ ) = const < 1. Assuming (3j^(r) 
(isothermally) without loss of generality, we obtain, 



(r)/c = P 3 and T JT (r) = T 3 
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Jo 

which can be applied to the cold plasma (electron- ion) jets for / P <1, electron-positron- ion jets for f p <^ 1, 
and electron-positron jets for f p ~ 1. 



3. ANALYTICAL SOLUTIONS: JET COLLIMATION AND TOROIDAL MAGNETIC 

FIELD 



The master equation (4) can be self-consistently solved for the total electron density n e -(r). In the 
case of [fbfib + (1 — f b )j3 r ] 2 > 1 — fc, implying that the confinement force of the self-generated magnetic field 
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exceeds the repulsive force due to the space charge, we find the solution in the form of the newly modified 
Bennett equilibrium: n e _(r) = n e _(r = 0)/[l + (r/r.j) 2 ] 2 . Here rj characterizes the radius of a filament of 
jet, which is given by 



rj = J Mj+MI f5) 

\j(i-f P ) 2 {[hh-(i-h)\M 2 -^-fc)Y 1 ! 

where \d = {[/fcTfc + (1 — /fc)T' r ]/47rn e „e 2 } 1 / 2 stands for the effective Debye sheath of electrons. For 
f c ~ 1 and f p ~ 1, the effective diameter of the charge- neutralized filament is approximated by 
dj = 2rj ~ <5/3- 1 (5/- 1 (4T , /n e _e 2 ) 1 / 2 . This scales as 

where = (1 - fp)' 1 , d^ 1 = [/ 6/ 3 & - (1 - fb)\Pr\}~\ and f = / 6 f 6 + (1 - f b )f r are the pair production 
rate, the current-neutral rate, and the effective thermal spread, respectively. When the return currents 
flow outside the filaments, i.e., their envelope and/or ambient medium, the current-neutral rate reduces to 

In Figure 1 for f c = 1, we show that the diameter of a filament swells as the pair production rate 
and return current increase and as the average density of the filament decreases. Physically, the increase 
of positron density leads to the seeming decrease of electron density, because only a small portion of 
electrons can take part in screening ions effectively, viz., 5f p n e - = (Z*)rii/ f c « {Z*)rn, and this results 
in significantly stretching the sheath radius. Moreover, the return current inside the filaments enforces 
to screen the magnetic fields, so that the filaments with the lower energy density of the magnetic fields 
prefer to be spread, to arrange the required magnetic confinement force against the thermal expansion and 
electrostatic repulsion. 

According to self-consistent analysis of particle orbits, the larger radius does not link with the 
larger current capacity (Honda 2000). The net current inside the filament sheaths will be strongly 
limited because of the orbital migration of electrons embedded in self-generated magnetic fields. That is, 
ij ~ 4eSf3cSf p n e ^rj £ 1.65/? 6 f h m e c 3 /e = 28.2/3 6 f b kA = i , where f b = (1 - ^)" 1/2 . This allows the 
parameter region of 1 > 5/35 f p J> 2.4T/f bm e c 2 . We conjecture that ejecta with the huge current /,/ > io 
(Appl & Camenzind 1992) necessarily split into many filaments, each carrying one "unit" current of ~ io, 
so that the number of the filaments is of the order of magnitude of Nf ~ Ij /io (Honda et al. 2000b). 
Note that the forward currents must be almost compensated with the backward return currents. The 
morphology seems to be consistent with filamentary structure including a backward flow discovered by the 
Highly Advanced Laboratory for Communications and Astronomy VLBI Space Observatory Programme 
survey (e.g., NGC 1275/3C 84; Asada et al. 2000), still one needs further study. 

Substituting the diffuse density profile into equation (3), we obtain the radial profile of magnetic field 
\Bg(r)\ = 2ire5j35f p n e -(r = 0)r/[l + (r/rj) 2 ]. At r — rj, the field takes the maximum value, to give for 
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Concerning the slowly decaying property of \Bg(r)\ cx r -1 , the actual diameter of a jet (i.e., a bundle 
of the numerous filaments) can be envisaged as Dj <~ 4(A r y/7r) 1/ ' 2 (_B| max 2 /87re a ) 1 / 2 (ij 3> dj, where e a 
is the energy density of the ambient medium. In contrast with the diameter Dj, the maximum field 
estimated in equation (7) does not depend on 5/3 and Nf. It is instructive to notice the relation of 
n e -(T b - l)m e c 2 > J B| max 2 /87r ~ (2/7r)n e _T, namely (3 = SirpajB 2 ~ 0(1) (Tajima & Shibata 1997), 
consistent with equipartition of excess kinetic energy of the jet which serves as a free energy source. 

In Figure 2, for f c = 1, we plot the range of possible magnetic fields for a given temperature and 
density of jets. The maximum by equation (7) may be a theoretical restriction to the field strength of the 
astrophysical jets, and close to the value required from the synchrotron self-Compton (SSC) model, which 
explains a double-humped appearance on spectral energy distributions (e.g., Mrk 501; Kataoka et al. 1999). 
On the other hand, from the spectral fitting by the modified synchrotron proton blazar (SPB) model, we 
expect a somewhat larger value (Miicke & Protheroe 2001), which is required from the strong synchrotron 
radiation, relating to stochastic acceleration of ultra-high energy cosmic rays. Quasi-perpendicular shocks, 
whose configurations are favorably set up in accordance with the present mechanism, can probably 
accelerate particles up to 10 — 100 EcV (10 19 — 10 20 cV). The parameter regions of other AGN jets 
(e.g., M87, 3C 273, Cyg A), which are still being argued at the moment, are also shown in Figure 2. More 
detailed and comprehensive observations are required in order to clarify the detailed activities in the tips 
and main bodies of jets (see e.g. Zavala & Taylor 2002). 



4. DISCUSSION AND CONCLUSIONS 

Finally, we remark that the present model is tolerant for the synchrotron cooling and the 
beam-plasma-type instability, which could be the major competitive processes. 

1. Let us suppose that beam electrons with a terminal energy of 0.1 — 1 GeV launched from a central 
engine lose their energies, by the synchrotron radiation due to their own magnetic fields, to 1 — 10 MeV, 
comparable to the transverse thermal spread. In this case, the synchrotron cooling time is of the order of 
T S yn ~ (10 2 — 10 4 )(0.1 gauss/£?) 2 yr, which corresponds to the propagation distance up to ~ 10 kpc for 
n e=F ~ 10 4 cm~ 3 and about a megaparsec for n eT <J 10 2 cm~ 3 . In addition, reacceleration processes of the 
electrons, if they work, further lengthen the propagation distance. 

2. Mutual coupling between beam and return currents as well as dissipation processes can cause various 
beam-plasma instabilities. Below we briefly explain the possible mechanism of the growth rate reduction in 
terms of the resistive hose instability as an example. For the square radial profile, single frequency seen by a 
beam particle can be resonant with the surface perturbation when = uj — kcj3b — (J3b5(35f p /T f,) 1 / 2 (D pe , where 
ij pe = (47rn e _e 2 /m e ) 1 / 2 denotes the mean plasma frequency. However, for the diffuse profile presented above, 
the growth rate is bounded, and there is no such single frequency f2 for which the entire beam is resonant with 
the wave (Uhm & Lampe 1980). According to the self-consistent Vlasov-Maxwell analysis, the "off-resonant" 
dispersion yields the most dominant wavenumber in the complex form of fc* = Tm(fci) ± Re(fc r ). We 
finally get the growth distance of Li = 27r/Im(fci) ~ TrcTd(3bS(3/[0.7fb(3b + (1 — /b)|/3 r |], where ~ irr 2 j<j/c 2 
stands for the magnetic diffusion time and a is the electrical conductivity, and the oscillation wavelength of 
A r = 27r/Re(fc r ) ~ (15c/w pe )[f bPb/(SpSf p )} 1,>2 . Note the relation of Li ^> A r , in contrast to Li ~ A r for the 
resonant case. When assuming the relativistic Spitzer conductivity a ~ 10 12 (T/10 10 i ; i') s _1 for T ^> m e c 2 
(Braams & Karney 1989), the distance Li can be expressed as 
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22 /icr 3 \ 2 ict 1 io~ 3 / t \ 2 io- 3 cm- 3 

for (3b — 1 and Sf p <C 1- The pairs expand the radius rj, to significantly increase the phase lag Td- 
Evidently, the screening effects by the pairs lower the growth rate. We also note that if there exist poloidal 
(longitudinal) magnetic fields superposed on the transverse fields (Gabuzda & Gomez 2001), the shear 
structure of the helical fields could be well established (Miyamoto 1989). They will play significant roles in 
stabilizing and guiding the jets (Davidson 1990) and be able to further stretch the distance L^. 

In conclusion, we have newly developed a generic model for collimation and transport of relativistic 
electron-positron-ion jets. A jet as a bundle of many filaments can be sustained owing to toroidal magnetic 
fields self-generated by the negatively charged stream. The magnetic field pressure balances with the gas 
pressure of the jet. The pair-screening effects and return currents expand the filament widths, and such 
expansion and the diffuse envelope of the filaments lead to strong suppression of the instabilities. In order 
to fully demonstrate observational results, three-dimensional relativistic Vlasov-Maxwell or electromagnetic 
particle-in-cell simulations with the larger spatiotemporal scale should be promoted in future. 

We acknowledge useful discussions with M. Kusunose, Y. Sentoku, and A. Mizuta. M.H. thanks the 
RIST and APR-JAERI for their hospitalities. This work was supported in part by the Grants-in-Aid of 
ITBL-Japan. 
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Fig. 1. — Effective diameter of a charge-neutralized electron-positron-ion filament dj vs. the pair production 
rate Sfp 1 = n e -/(n e - — n e+ ) {top) and the current-neutral rate S/3' 1 {bottom)-, for an effective temperature 
of T = 10 10 K and average electron densities of n e _ = 10~ 3 , 1, and 10 3 cm -3 . We have chosen the typical 
parameters of the current-neutral rate and the pair production rate: <5/3 _1 = 1 {top) and Sfp 1 = 10 3 {bottom), 
respectively. For further explanation see text. 
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Fig. 2. — Maximum possible magnetic fields Bg max for a given electron density n e _ and effective temperature 
(Tio = T/10 10 K) of jets. The solid and dotted lines show the maximum fields for the electron-ion jets 
(fp = 0) an d the electron-positron jets (f p — 1), respectively. The shaded areas indicate the allowable 
parameter regions of some well-known AGN jets: Cygnus A (Carilli et al. 1998), Virgo A/M87 radio lobe 
(Owen, Eilek, & Kassim 2000), 3C 273 (Aharonian 2001), and Markarian 501 (predicted by the SSC model 
[Kataoka et al. 1999] and by the modified SPB model [Mucke & Protheroe 2001]). Note that all quantities 
are in the ion rest (jet) frame. 



